Preferential attachment growth model and nonextensive statistical mechanics 
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We introduce a two-dimensional growth model where every new site is located, at a distance r from 



the barycenter of the pre-existing graph, according to the probability law 1/r 



{ac > 0), and 



is attached to (only) one pre-existing site with a probability oc ki/r"^ {ua > 0; ki is the number 
of links of the i*'' site of the pre-existing graph, and its distance to the new site). Then we 
numerically determine that the probability distribution for a site to have k links is asymptotically 
given, for all values of ac, by P{k) oc where = [1 -f (1 — g)x]^''''^~''' is the function 

naturally emerging within nonextensive statistical mechanics. The entropic index is numerically 
given (at least for oa not too large) by g = 1 -f- (l/3)e~°'^^^ and the characteristic number 
of links by K ~ 0.1 -I- 0.08 qa. The a a ~ particular case belongs to the same universality class 
to which the Barabasi- Albert model belongs. In addition to this, we have numerically studied the 
rate at which the average number of links {ki) increases with the scaled time t/i; asymptotically, 
{ki) oc {t/i)^ , the exponent being close to j3 = ^(1 — ua) for < < 1, and zero otherwise. 
The present results reinforce the conjecture that the microscopic dynamics of nonextensive systems 
typically build (for instance, in Gibbs F-space for Hamiltonian systems) a scale-free network. 



Among the subjects that are being studied intensively 
nowadays in statistical physics, there are two, namely 
nonextensive statistical mechanics (see Q for a review) 
and networks ,2i], in particular scale-free networks 0, 
which receive special attention in connection with com- 
plex systems 0,|EIEI3- Could these two topics be inti- 
mately related? This would not be so surprising after all, 
given the fact that both research lines frequently address 
similar types of natural and artificial systems, in physics, 
economics, chemistry, biology, linguistics, social sciences 
and others. In fact, such a connection has already been 
conjectured in several occasions, e.g., Q (Preface and 
Chapter 1) and j8]. In the present paper we propose a 
growth model, on which we exhibit and quantitatively 
explore this connection. 

Let us consider a continuous plane. We shall construct 
a single connected network of sites (or nodes or vertices) 
and links (or bonds or edges) by gradually (sequentially) 
making it grow. We first fix one site {i = 1) at some arbi- 
trary origin of the plane. The second site {i = 2) is ran- 
domly and isotropically chosen at a distance r distributed 
according to the probability law Pq (r) c>c 1/ r^+"<5 [ac > 
0; G stands for growth). This second site is then linked 
to the first one. To locate the third site {i — 3) we move 
the origin to the barycenter of the two first sites, and 
apply again the distribution Pair) from this new origin. 
This third site is now going to be linked to only one of 
the pre-existing two sites. To do this, we use an attache- 
ment probability pA oc ki/r"^ {aA > 0; A stands for 
attachment), where is the distance of the newly ar- 
rived site to the i*'* site of the pre-existing cluster, and 
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FIG. 1: Typical iV = 250 network for (qcQa) = (1,0). The 
starting site at {X, Y) — (0, 0) is indicated with a larger circle. 
Notice the spontaneous emergence of hubs. 



the connectivity ki is the number of links already arriving 
to the same i*^ site (at the present stage, ki — k2 — 1). 
The growth-attachment process is sequentially repeated 
as long as we want; the sites have unit radius, and no 
new site is admitted which overlaps with a previous one. 
If we denote with N the total number of sites of the clus- 
ter, it immediately follows that the linking of the newly 
arrived site (i — N) is done with the probability 
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It is clear that the dynamics of this model makes the 
arriving sites to have preferential attachment to the pre- 
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FIG. 2: Connectivity distribution for qa = 1 and typical 
values of ac; 2000 realizations of iV = 10 000 networks . 



vious sites that already have many links (hubs), as long as 
they are not too far. This competition (already explored 
in Q for the particular case of uniform distribution of 
sites within some limited region) between connectivity 
and (Euclidean) proximity is less pronounced when a a is 
close to zero, and completely disappears only at = 0. 
For this particular case, one expects behaviors consistent 
with the Barabasi- Albert model 0, which has topology 
but no metrics. In the present paper, we focus on two 
main aspects: (i) the stationary-state connectivity distri- 
bution P{k) associated with the number of sites that have 
k links in the N oo limit; (ii) the time dependence of 
the average number (ki) of links, more precisely how (ki) 
grows with the scaled time t/i (t > i), particularly in the 
limit t/i ^ oo (see 0)- 

Typical networks obtained with this model are indi- 
cated in Fig. 1. Our numerical results for P{k) are in- 
dicated in Figs. 2 and 3. We illustrate in Fig. 2 the 
fact that P{k) does not depend from ac for any given 
value of a A- The exponent ac controls the metrics of the 
emerging cluster, but has no influence on the connectiv- 
ity distribution. In contrast, this distribution is greatly 
influenced by the exponent aA, as illustrated in Fig. 3. 
We can verify that all our examples are very well fitted 
with the form 

P(fc) = P(0)e-'=/% (2) 
where the q-exponential function is defined as follows 0, 
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K > is a characteristic number of links. In Fig. 4 we 
show the aA dependences of q and k. In particular, q{aA) 
does not exhibit a critical value of ua above which a 
different regime could emerge. However, unless a detailed 
analysis (out of the aim of the present letter) is done of 
the finite-size effects, such a possibility should not be 
excluded. 



FIG. 3: Connectivity distribution for typical values of oa (we 
have used ac ~ 2 but we recall that this value is irrelevant). 
Points are our computer simulation results; continuous lines 
are the best fits with the q-exponential function indicated 
in Eq. (2). (a) log — log representation; (b) ln4/3 —linear 

representation, with In^ a; = ^ i-g^ ■ (^) —linear repre- 
sentation, where, for each value of aA, we have used its corre- 
sponding value of q. We have used 3 different representations 
to improve comprehension. 
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FIG. 4: Values of q and k used in the best fits indicated in Fig. 
3. The solid curves are: (a) g = 1 + (l/3)e-°-^^'' (V ao); 
(b) K ~ 0.083 + 0.092 aA (V aa)- 

We must now address a relevant point: Why have we 
fitted our curves with (/-exponentials? The reason lies on 
the conjectured connection with nonextensive statistical 
mechanics [13 (see 0, [lll | for reviews), a theory which 
we briefly illustrate now. Let us consider the following 
generalization of the Boltzmann-Gibbs entropy Sbg = 
-/ dkP{k)\nP{k): 

Sq = (q eU; Si= Sbg) (4) 

9-1 

If wc optimize this entropy with the constraints 
JdkP{k) = 1 and {J dkk[P{k)]'i}/{J dk[P{k)]''} = K 
we obtain straightforwardly 12] 

-r3rAk-K) 

m = , ' . , (5) 
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FIG. 5: (a) Time dependence {t = A'') of the average number 
(over 2000 realizations) of links for typical values of ua for 
sites i = 1, i = 5 and i = 95 (compare with Fig. 2(c) of '^). 
We have used aa = 2; (b) aA-dependance of /3 (the straight 
line (3 — ^{1 ~ a a) could be the exact answer). 

where (3q is simply related with a Lagrange parameter. 
This expression precisely coincides with Eq. (2) through 
the identification 

Before proceeding with our numerical results, let us 
briefly mention what indications make us to beUeve that 
a close connection, and not just a mere functional co- 
incidence, might exist between the present preferential- 
attachement growth model and the thermostatistical sys- 
tems addressed by the entropy (4). First, a growth model 
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involving, like the present one, preferential attachment 
has been proposed and both analytically and numerically 
discussed in . In this model the connectivity distribu- 
tion is analytically shown to be precisely of the form (2) 
(although written in a slightly different manner). The 
corresponding entropic index is given Q (Chapter 1) by 
q = [2m(2 - r) + I - p - r]/[m(3 - 2r) + 1 - p - r], 
where {rn,p,r) are parameters of the model. Second, 
Lennard- Jones small clusters (with up to 14 atoms) have 
been numerically studied recently [l4| . The distributions 
of the number of local minima of the potential energy 
with k neighboring saddle-points in the configurational 
phase space can, although not mentioned in the original 
paper , be quite well fitted with g-exponentials with 
g = 2 Third, the present model generates structures 
that are scale-free (see, for instance. Fig. 1). Consis- 
tently, there is plethoric evidence of the connections of 
the nonextensive entropy Sq, and of its associated statis- 
tics, with hierarchical and (multi)fractal structures (see, 
for instance, |0, Ell)- It is clear that none of 

these three features constitutes a proof; however, the set 
of them suggests the quite plausible scenario that the 
present growth model basically satisfies the dynamical 
requirements for the nonextensive concepts to be appli- 
cable. 

Let us now present our results concerning the rate at 
which the number of links of a given site increases with 
time. A basic quantity is the average (over 2000 realiza- 
tions) number (ki) of links at a given time, more precisely 
as a function of the scaled time t/i j5|. Our results are 
illustrated in Fig. 5. We numerically find that, for all 
values of (0^,00), (h) oc {t/i)'^'^'^^\ with /?(0) = 1/2. 



Let us summarize the present paper. We have intro- 
duced a growth model which has both preferential at- 
tachment and metrics. Every new site of the cluster has 
to "decide" to which one of the pre-existing sites will link 
itself. In this stochastic choice, there might be compe- 
tition between the "popular" sites (those that already 
have many links) and the nearby sites. This competition 
is more pronounced for increasing value of the exponent 
a A', it disppears for a a = 0. The connectivity distribu- 
tion and the rate of increase of links with time are sub- 
stantially influenced by a a- The a a = model belongs 
to the same universality class to which the well known 
Barabasi-Albert model belongs. In addition to these re- 
sults, we have shown that the connectivity distribution 
is (numerically) given by the g-exponential function that 
emerges naturally in the frame of nonextensive statistical 
mechanics. This fact supports the conjecture that the 
typical occupation of the accessible phase space of many 
nonextensive dynamical systems might be, in relevant 
stationary states, scale-free. This is in notorious variance 
with the occupation of standard, extensive, isolated dy- 
namical systems, which tends to be uniform, and whose 
equilibrium state is that prescribed within Boltzmann- 
Gibbs statistical mechanics. 
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